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Abstract 

We discuss effect of nuclear shadowing in neutrino deep-inelastic scattering in terms of non 
| perturbative parton model. We found that for small Bjorken x and large Q 2 the structure 

function F$ is shadowed in nuclei about two times as stronger as F^. The underlying 
reason and phenomenological aspects of this observation are discussed. 
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m ! 1 Introduction 

The small x physics is an important topic of high energy studies nowadays. We note 
. in this respect that data is usually collected for nuclear targets and it is important to 
. consider nuclear effects in deep inelastic scattering (DIS). The systematic reduction of the 
CN \ nuclear structure functions FA/A with respect to the nucleon F% has been established 
experimentally at CERN Q and Fermilab || in the region of small Bjorken x, so called 
Q\ • nuclear shadowing phenomenon. This phenomenon has been extensively discussed in the 
c-| ■ literature within different theoretical approaches (for a review see ,e.g., ref.0). 
Ph. Apart from F2, it is important of course to study the nuclear effects in the other 

Ph! observed structure functions. An interesting object in this respect is the structure function 
F 3 which is measured in (anti)neutrino DIS. We recall that at large momentum transfer 
F 3 describes the valence quark distribution in the target. The main concern of the present 
paper is to study the nuclear effects in neutrino DIS in the small x region and compare 
the shadowing effect for the structure functions F2 and F$, i.e. for the sea and the valence 
quark distributions. 

For the discussion of nuclear physics at DIS it is convenient to choose the target rest 
frame. To clarify the space-time picture of DIS in this reference frame it is convenient to 
argue using the time-ordered ('old') perturbation theory. We recall that DIS cross section 
is determined by the imaginary part of the Compton scattering amplitude of the virtual 
intermediate boson from the target. Two typical time-ordered diagrams which contribute 
to the Compton nucleon scattering amplitude are drawn in Fig.l. 
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a FIG.l b 

Time-ordered diagrams (t\ < t2) for the Compton scattering. 



1 Based on the talk given at XIV International Seminar 'Relativistic Nuclear Physics And Quantum 
Chromodynamics ' Dubna, August 17-22, 1998. To be published in proceedings. 
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Fig. la describes the scattering process which goes via absorption and re-emission of the 
virtual boson by a quark bound in the target. Fig. lb shows the process where the virtual 
boson first decays into a gg-pair which then interacts with the target. By comparing the 
energy denominators of these two diagrams one can observe that at small x the diagram 
in Fig. lb dominates over the one in Fig. la. This observation leads to the qualitative 
space-time picture of the Compton scattering at small x which consists from two time- 
ordered processes. First, the virtual boson splits into a qq (or hadronic) fluctuations 
before arrival at the target. Then these states propagate and interact with the target. 
A typical propagation length of a qq fluctuation can be estimated from the uncertainty 
principle and at high Q 2 is of order L ~ (2Mx) _1 with M the nucleon massQ M. For 
2; < 0.1 the qq fluctuation lives long enough to develop a chain of multiple interactions 
with bound nucleons. This will lead to nuclear shadowing effect which we discuss in detail 
in Sec.^|. Before we come to the discussion of nuclear modifications in the DIS, in Sec.|2| 
we set up the model for the DIS from the isolated nucleon and fix some parameters of the 
model by comparing to data. 

2 DIS in non perturbative parton model 

We consider the neutrino and antineutrino charged-current DIS from the nucleon. It is 
supposed that Q 2 is sufficiently high to neglect the target mass effects. We neglect also 
possible final state interaction effects and focus on the leading twist contribution to the 
DIS which is given by the handbag diagram of Fig. 2 (note that this diagram is treated in 
terms of a covariant approach which includes of course all time-ordered processes shown 
in Fig.l). 
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FIG.2 

Compton scattering amplitude from the nucleon to the leading order in Q 2 . 

In the axial gauge the Compton amplitude for the scattering of virtual IV-boson from the 
target reads (to be specific we consider the Compton scattering of W /+ -boson): 
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where q is momentum transfer and Xa{k, P) is the propagator of the quark fields of kind 
a in the target with momentum P, 



Xa(k,P) 



d ^e* k <(P\TU a (OM0))\P)- 



(2) 



L is also a typical time between absorption and emission of the virtual boson in Fig.l 
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The dots in Eq.(|T]) denote contributions from the s- and c-quarks which repeat the first 
two terms with the substitutions u — > c and d — > s. 

The structure functions are given by the imaginary part of the Compton amplitude 
and can be found from Eq.fll]) by applying appropriate projection operators. In the leading 
order in 1/Q 2 we find that the structure functions F 2 and F 3 are given in terms of the 
quark light-cone distributions by the standard parton model equations (neglecting a small 
effect due to Cabibo mixing angle), 

F 2 v {x) = 2x(d(x) +u(x) + s(x) + c(x)), (3) 
Fg(x) = 2 (d(x) - u{x) + s(x) - c(x)) . (4) 

where the parton distributions are expressed in terms of the quark propagator as follows, 
JaK ' J (2tt) 4 2P-q \ P-qJ ' V ; 

Qa(x) = f a (x), q a (x) = -f a (-x). (6) 

Note that the quark and antiquark distributions come from the direct and the crossed 
terms of the Compton amplitude respectively. It follows from Eq.fll]) that the structure 
function F\ is not independent; it is given by the Callan-Gross relation, F 2 (x) = 2x Fi(x). 

In QCD scaling is violated and the parton distributions depend also on Q 2 . The Q 2 
dependence comes through the dependence of the quark correlator on the normalization 
point fi 2 which is taken to be Q 2 in our discussion. 

For our purposes it will be convenient to write the parton distributions as the integrals 
over the invariant mass spectrum of the spectator states ||. To this end we assume 
the analytic properties of the quark correlator (0) to be similar to the ones of hadronic 
amplitudes, i.e. it is an analytic function of the invariant variables s = (P — k) 2 , u = 
(P + k) 2 and k 2 . For real s and u the quark amplitude has a right-hand cut in the variable 
s, a left-hand cut in the variable u and singularities at k 2 > 0. In order to make use of 
these analytical properties of in the loop integral (|5]), it is convenient to parameterize the 
loop momentum k in terms of the light-cone or the Sudakov variables. Acting along these 
lines one finds that the distribution function f(x) vanishes outside the physical interval 
\x\ < 1. For < x < 1 the distribution function f(x) (or the quark distribution q(x)) 
is given by a dispersion integral along the cut in the s-channel, while for — 1 < x < 
the ti-channel cut is relevant. We separate from the quark correlator (0) external quark 
propagators and write the result in terms of the quark-nucleon, T q ^, and the antiquark- 
nucleon, Tg N , scattering amplitudes. In terms of these amplitudes the small x part of the 
quark and antiquark parton distributions read as follows (to simplify notations we drop 
here explicit dependence on the quark kind a) 

q{x) = (2^J dS J dk (k*-m*)* { ~ X + - - (7) 

<iu) = 7-—,l<i«l de lv ff {u ':: ■ ( - > + , i, (si 




N r . r , ^maxK^) ImToArfw, k* 

v -. v . . v .. .4) 

Here A^ c is the number of colors, m q is the quark mass. The integrations run over the 
spectrum of intermediate states and over the quark four-momentum squared k 2 with 
k 2 asx (s, x) = x[M 2 — s/(l — x)) being the kinematical maximum of k 2 for the given s 
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and x. We note that the quark distribution q(x) is determined by the giV-scattering 
amplitude, in accordance with intuitive discussion based on time-ordered diagrams of 
Fig.l. For the antiquark distribution q(x) the giV-amplitude is relevant. The reader can 
find the derivation of Eqs. (|7]||) in ref. 0. 

Eqs.([7|||) establish the correspondence between the x-dependence of parton distribu- 
tions and the region of invariant masses of the residual quark system. The basic assump- 
tion of the non-perturbative parton model is that the parton amplitudes vanish at large 
virtualities so that the integrals in Eqs. (|7|,|D are finite and saturated in the region of fi- 
nite k 2 . At small x the k 2 is finite even for large s ~ M 2 /x and therefore the parton 
distributions are sensitive to the large-s part of the spectral densities. 

In ref. (see also 0) we have considered a simple model for the spectrum of the 
residual quark system which incorporates two scattering mechanisms shown in Fig.l. It 
was assumed that the spectrum of spectator states is divided into low-mass and high-mass 
parts sharply at some separation parameter sq. The physical origin of the low- mass part 
of the spectrum is due to the Compton scattering of the intermediate boson from bound 
quarks in the target, Fig. la. This process is of major importance at intermediate and 
large x. The high-mass part of the spectrum (s > sq) gives rise to the low-x part of the 
structure functions. In this region the spectrum is dominated by the process shown in 
Fig.lb. 

In order to calculate the quark and anti-quark distributions one has to specify the 
amplitudes T q ^ and T q N- We follow here ref. and extract the quark-nucleon amplitudes 
from the observed high-energy hadron scattering amplitudes. The idea is to approximate 
the amplitudes T q n and T q N by the constituent quark-target scattering amplitudes which 
we obtain from the additive quark model of high-energy scattering. We are concerned 
here with the scattering from the isoscalar target. The isoscalar parts of the (anti)proton- 
nucleon forward scattering amplitudes, T pN = (T pp + T pn )/2 and T pN = (T pp + T pn )/2, can 
be written in terms of the (anti) quark-nucleon amplitudes as, | 

T pN (S)/S = (T qN (yS,k 2 )/yS), 

T pN (S)/S = (T qN (yS,k 2 )/yS). (9) 

Here S is the nucleon-nucleon invariant mass squared, k is the four-momentum of a 
constituent quark in the target and y is a fraction of the target light-cone momentum 
carried by the constituent quark. The averaging is taken over the quark distributions in 
the proton which is normalized to the number of quarks. 

An approximate solution of Eqs.@ can be obtained by evaluating the qN and qN 
amplitudes at averaged values of y and k 2 for constituent quarks in the proton, y ~ 0.2, 
k 2 ~ —0.3 GeV 2 . The s-dependence of the (anti)quark amplitudes is then fixed by the 
energy dependence of the observed pN and pN scattering amplitudes. In the Regge 
approach these amplitudes can be approximated by the exchange of the Pomeron P and 
two Regge poles corresponding to the scalar (/) and the vector meson (u) trajectories [Q, 
T p n = P + f — oj, TpN = P + f + uj. Eqs. (|D thus determine the quark-nucleon Regge 
poles residues in terms of the nucleon-nucleon ones. The exchange of the Regge trajectory 
with the intercept cxr leads to the imaginary parts of the amplitudes rising as s aR . The 

3 In fact if the isospin is the exact symmetry one finds T u n = T c in — T q N and T u n = T^ N = T q N- 
We have used this observation implicitly by writing Eq.(||) in terms of only two isoscalar quark-nucleon 
amplitudes. 
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Pomeron exchange dominates both q(x) and q(x) at small x which rise asymptotically as 
x~ ap , as one can see from Eqs.(|7|||). The Pomeron and the scalar /-Regge pole cancel 
out in the difference q — q (the valence quark distribution). Note however that the valence 
quark distribution remains finite at small x with the asymptotics determined by the u 
Regge trajectory, q(x) — q(x) ~ x~ a " . In the numerical estimates we use the parameters 
of best fit of ref.|| with the Regge poles intercepts ap = 1.08 and af = a u = 0.55. 

In the model of ref.|| it is assumed that the coupling of the Regge poles to quarks 
is independent of its kind, and thus the different parton distributions are described by 
two amplitudes T qN and Tq N ^ In order to have the integrals in Eqs.(|7]|D convergent 
the amplitudes should fall sufficiently fast for large quark virtualities k 2 . We assume the 
factorized s- and fc 2 -dependencies for every Regge pole in the quark-nucleon amplitude 
which we parameterize in a simple form (1 — k 2 /Ap)~ nR . The cut-off and exponent 
parameters as well as the scale parameter sq of the spectrum of spectator states are 
chosen to fit the charged-lepton and neutrino data for the structure functions F2 and xF% 
at average momentum transfer Q 2 between 5 and 10 GeV 2 (for more details see ref.|J). 

A comment on the scaling violation is in order. In our approach we fix the parameters 
of the spectrum of spectator states to describe the measured structure functions at some 
fixed Q 2 and therefore effectively incorporate the Q 2 evolution effect in the phenomeno- 
logical parameters of the spectrum. A more consistent approach which accounts for the 
scaling violation effect would be to choose a low resolution scale, Q 2 ~ 1 GeV 2 , and fix 
the parameters of the spectrum at this scale. Then at larger Q 2 the parton distributions 
are obtained by applying the Q 2 evolution equations. This program in fact has been 
attempted in ref.]7|]. 

3 Nuclear Shadowing 

We now turn to discussion of DIS on nuclear targets in the small x region. The physical 
origin of nuclear corrections in the small x region is the coherent interaction the propa- 
gating qq pair with bound nucleons. In the parton model regime, which is discussed in 
this paper, the qq pair is highly asymmetric in the target rest frame, i.e. one of the pair 
constituents carries practically all the transferred momentum. It is commonly assumed 
that the fast quark only weakly interacts with the target and begins to radiate soft gluons 
and hadronize well after the nucleus. The other participant, the wee (anti) quark, has a 
finite momentum and interacts in the nucleus with normal hadronic cross-section. 

It is instructive to compare different time and length scales involved in this process. 
The onset point of coherent nuclear effects x Q is usually estimated by the comparison of 
the propagation length L ~ (2Mx) _1 of the qq pair with the averaged distance between 
bound nucleons in the nucleus. If we take for estimates the averaged nuclear density 
p = 0.17 fm -3 in the central nuclear region then x Q ~ 0.1. We note however that this 
is only a kinematical condition which allows for coherent interactions of the propagating 
qq pair with several nucleons for x < x a . It is quite clear that the region of developed 
nuclear shadowing should depend on (anti) quark-nucleon interaction. Indeed if a is the 
averaged (anti) quark-nucleon total cross section then the (anti) quark mean free path 
in the nucleus I = (per) -1 . The multiple scattering effect becomes significant when the 

4 We note in this respect that the parton distributions depend on their masses, as one can see from 
Eqs.(fnJq). Therefore the contribution from heavy quarks is suppressed by their large masses. 
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coherence length L exceeds the (anti)quark mean free path in the nucleus I. This happens 
when x < pa/2M « 0.02, where we have taken for estimates a = (7^/3 ~ 13 mb. In this 
case the main part of the incoming flux of qq pairs is absorbed at the nuclear surface in 
the layer with the thickness of order I and thereby the inner nucleons are screened. 

In order to quantify the discussion we consider the quark and antiquark nuclear scat- 
tering amplitudes in forward direction, T qA and T qA , and apply the Glauber-Gribov mul- 



tiple scattering expansion [ ID , O ] assuming that formalism can be taken over to off-shell 
interaction, 

T A = AT N + 5 {2) T A + ... (10) 
i5^T A (s } k 2 ) fiT N (s,k 2 ,ty ~ 



2s 



A(A-l) [ N[S 2 [ s - M jd 2 bdzdz'e^-*')pM(b,z,z'). 



z'>z 



Here we have written explicitly only the first and the second terms in the expansion. The 
first term accounts for the incoherent scattering from bound nucleons while the second 
one is a correction due to the double scattering. The dots denote higher order scattering 
terms. The notations are as follows: k = (ko, 0j_, k z ) is the momentum of incoming off- 
shell (anti)quark, q\\ = (m 2 — k 2 )/2\k z \ is the longitudinal momentum transfer in the 
scattering of the (anti) quark from the nucleon, t = —q 2 , s is the (anti) quark- nucleon 
center-of-mass energy squared and p( 2 >(b,z,z') is the two-nucleon density matrix. In 
our estimates we neglect effects due to inelastic excitations of the projectile particle in 
intermediate states as well as the effects of nucleon-nucleon correlations in the nucleus 
and assume that p^ factorizes into two one-particle densities normalized to the unity, 
p( 2 '(b, z, z') = p(b, z)p(b, z'). In the numerical evaluations we use the Gaussian shape 
for the densities p(b, z) with the scale parameter fixed by the nuclear root-mean-square 
radius. 

We use Eq.(|l(]) to calculate (anti) quark- nuclear scattering amplitudes and then eval- 
uate nuclear corrections to the quark q(x) and antiquark q(x) parton distributions in the 
isoscalar nucleon by Eqs.([7|||). We keep the terms up to the 4th-order in the multiple 
scattering series (0). In Fig.3 we plot the ratios R 2 = F^/AF^ and R 3 = FA/AFf for 
the calcium nucleus for the structure functions averaged over neutrino and antineutrino 
fluxes, F 2 = (F 2 u + F%)/2 and F 3 = (F£ + Fg)/2 (the structure functions are evaluated for 
the isoscalar nucleon). The double scattering correction is negative^ and brings the main 
contribution to the multiple scattering series. However the share of the 3d-order and the 
4th-order terms is also sizable, about 25% of the double scattering term in F 3 and about 
15% in the case of F 2 for the calcium nucleus. The sum of the 3d-order and the 4th-order 
terms is positive and partially compensates the double scattering term. Note also that 
the share of higher order terms increases in heavy nuclei. 

One can observe from Fig.3 that the nuclear shadowing effect for F 3 (i.e. for the 
valence quarks) is about 2 times stronger as the one for F 2 (i.e. for the sea quarks). To 
clarify the underlying reason for this factor of 2 we recall that the structure functions F 2 
and F3 are determined by the sum and by the difference of the antiquark and the quark 
scattering amplitudes respectively. At high energy the Pomeron pole dominates both the 
quark and the antiquark scattering amplitudes. The difference of the amplitudes, which 



5 This is well seen fromEq.© if we neglect the effect of g|| as well as a small real part of the scattering 
amplitudes. 
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is determined by the cu-Reggeon contribution, is small as compared to their sum. Since 
the double scattering correction in Eq. ([lO|) is bilinear in the quark amplitude, we observe 
immediately that the relative correction to the difference of the quark and the antiquark 
cross sections is two times as bigger as the corresponding correction to their sum.[] 




0.0001 



FIG.3 

Ratios i?2 and R3 calculated for the 40 Ca nucleus for small x and fixed Q 2 = 5 GeV 2 (see 
text). The data points for the ratio i?2 are from NMC [jl]] (filled boxes) and E665 [Q] (open 
circles) experiments with muon beams. 

We found a sizable correction to the Gross-Llewellyn-Smith sum rule for heavy nuclei 
due to nuclear shadowing. For the calcium and iron nuclei the corrections are 



0.3 



dx 



10- 



A 



-0.11 ( 40 Ca) 
-0.12 ( 56 Fe) 



It is usually believed that the GLS sum rule should not be renormalized by nuclear effects 
in the leading in 1/Q 2 order, because in this order the GLS integral counts the baryon 
number of the target. It is challenging therefore to look for a mechanism which would 
compensate a negative correction due to nuclear shadowing ('antishadowing'). Such an 
antishadowing may come from x > 0.1 region. We note in this respect that the Fermi- 
motion and nuclear binding corrections, which are of importance at large x, do not violate 
the GLS sum rule [I2|. A possible 'compensating correction' could be due to off-shell 



modification of the bound nucleon structure functions. 

In Fig. 3 for comparision we put also the data points for R2 for 40 Ca nucleus from NMC 
I] and E665 experiments with muon beams. We note that data is taken for different 
Q 2 for every rc-bin and, due to the nature of the fixed target experiments, small x is always 
correlated with small Q 2 (for events with x < 10~ 2 presented in Fig. 3, averaged Q 2 < 
1 GeV 2 ). In the small Q 2 region where scaling is violated, the applicability of the leading 
twist calculation is questionable. We therefore need a model which provides a smooth 
transition from scaling to nonscaling region. In 0] we have considered a model which 



The author is grateful to N. N. Nikolacv for useful discussion on this point. 
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interpolates the structure function Ftf (measured in muon induced reactions) between 
scaling and nonscaling regime where the latter was described in terms of vector meson 
dominance model [H| (VMD). It was found in the framework of this two "phase model" 
that the VMD part of F 2 is shadowed strongly than its scaling part. This lead to the 
enhancement of the overall nuclear shadowing effect. A similar observation has been done 



in a recent analysis of nuclear shadowing effect in neutrino induced DIS [[14] . In conclusion 
we notice the importance of a small Q 2 analysis for the structure function F 3 . We recall 
in this respect that F 3 describes the axial-vector-vector current transitions in the target 
under the charged current and the studies of this effect may bring new insights into the 
problem of chiral symmetry violation in nuclear environment. 

The author is grateful to W. Melnitchouk and N. N. Nikolaev for useful discussions and 

to the organizers of the Seminar for warm hospitality. 
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